arXiv:1505.00613v2 [hep-lat] 12 May 2015 


Evidence for the existence of udbb and the non-existence of ssbb and ccbb tetraquarks 

from lattice QCD 


^^^Pedro Bicudo, Krzysztof Cichy, ^^^Antje Peters, ^^^^Bjorn Wagenbach, and ^^^^Marc Wagner 

CFTP, Dep. Fisica, Instituto Superior Tecnico, 

Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 
Johann Wolfgang Goethe- Universitdt Frankfurt am Main, Lnstitut fiir Theoretische Physik, 
Max-von-Laue-Strafle 1, D-60438 Frankfurt am Main, Germany and 
Adam Mickiewicz University, Faculty of Physics, Umultowska 85, 61-614 Poznan, Poland 

We combine lattice QCD results for the potential of two static antiquarks in the presence of two 
quarks qq of finite mass and quark model techniques to study possibly existing qqbb tetraquarks. 
While there is strong indication for a bound four-quark state for qq = {ud — du)fV2, i.e. isospin 
7 = 0, we hnd clear evidence against the existence of corresponding tetraquarks with qq € {uu, (ud + 
du)j\/2, dd}, i.e. isospin I — 1, qq = ss and qq = cc. 

PACS numbers: 12.38.Gc, 13.75.Lb, 14.40.Rt, 14.65.Fy. 


I. INTRODUCTION 

Exotic hadrons have been proposed many years ago. 
As soon as quarks were found in the sixties, it became 
clear that systems more complex than standard mesons 
{qq states) and baryons {qqq states) could possibly exist. 
However, exotic hadrons are very elusive systems. Con¬ 
firming their existence or non-existence still remains one 
of the main challenges of particle physics. 

Frequently discussed exotic hadrons are tetraquarks 
Bia, which are four-quark bound states composed of two 
quarks and two antiquarks. There are several hadronic 
resonances which are tetraquark candidates. Among 
them are the light scalar mesons cr, k, /o(980) and 
ao(980) as well as the heavier mesons D*q and Dgi- How¬ 
ever, these systems have quantum numbers also consis¬ 
tent with a standard qq structure and their masses are 
not too different from what is expected in a qq picture. 
Thus, it is hard to rigorously argue that they are in¬ 
deed predominantly tetraquarks. On the other hand, 
there are also candidates which have quantum numbers 
or masses typical for tetraquarks, but not for standard qq 
mesons. For example 71^'^ has exotic quantum numbers 
= 1- ■+ or Zf and masses and decay products 
strongly suggest hidden cc or bb pairs, while their electri¬ 
cal charge ±1 indicates isospin 7=1. While the evidence 
for TTffP is not conclusive and the Z^ claimed by the 
BELLE collaboration [3] remains to be confirmed by dif¬ 
ferent experimental collaborations, the Z^ has received 
a series of experimental observations by the BELLE col¬ 
laboration mis], the Cleo-C collaboration [^, the BESHI 
collaboration [THII] and the LHCb collaboration m- 
Nevertheless, the Z^ would profit by more comprehen¬ 
sive measurements of its decay channels. We expect the 
existing and future experimental collaborations to con¬ 
tinue the study of present tetraquark candidates and to 
possibly also discover further ones. 

The theoretical study of tetraquarks is crucial to con¬ 
firm and correctly interpret corresponding experimen¬ 


tal observations and could as well provide information 
in which channels tetraquarks may be found. How¬ 
ever, tetraquark studies face a number of difficulties, e.g. 

(1) tetraquarks are usually open to meson-meson decay, 

(2) tetraquarks are complex relativistic four-body sys¬ 
tems, (3) quark models still fail to reproduce sectors of 
standard hadronic spectra and, thus, are not yet suffi¬ 
ciently well calibrated to reliably predict tetraquarks. 

In this work, we study the existence/non-existence of 
tetraquarks with two heavy bottom antiquarks bb. To 
this end, we use potentials of two static antiquarks in the 
presence of two quarks qq of finite mass, which we com¬ 
pute using lattice QCD. We extend recent studies of qqbb 
tetraquarks min] , where qq G {(ud—du)/-\/2 , uu, (ud-h 
du)/\/2, dd}, to similar systems with heavier quarks, 
qq = ss and qq = cc. In the future, we also plan to 
extend our investigations to the bb tetraquarks claimed 
by the BELLE Collaboration [3] . Such tetraquarks with 
a bb pair are, however, rather difficult to study with lat¬ 
tice QCD, since they couple to several decay channels. 

We avoid some of the technical difficulties of study¬ 
ing tetraquarks following a strategy already identified in 
the eighties [15]. We search for bound states rather than 
for resonances, to avoid open decay channels. Moreover, 
by using bb potentials obtained by lattice QCD com¬ 
putations, we largely avoid the calibration problem of 
quark models. Very heavy antiquarks such as b allow for 
the Born-Oppenheimer approximation |16j . For the two 
lighter quarks qq, the heavy antiquarks bb can be approx¬ 
imated as static color charges, which allows to determine 
the light quark energy using lattice QCD. On the other 
hand, once the energy of the light quarks qq is deter¬ 
mined, it can be utilized as an effective potential for the 
heavy antiquarks bb. 

Our lattice QCD computation goes beyond computa¬ 
tions with four static quarks, which show a clear evi¬ 
dence for four-body tetraquark potentials [m [iH] and 
tetraquark flux tubes [191 [20]. On the other hand, lattice 
QCD computations with four quarks of finite mass are 
extremely difficult and have found neither evidence for 
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charmed tetraquark bound states with udcc flavor |21j 
nor for resonances in the Zc family |22j . 


This paper is organized as follows. In section [ITj we 
briefly review the quark model and discuss qualitative 
expe ctat ions regarding qqbb four-quark systems. In sec¬ 
tion III we discuss the lattice QCD computation of bb 


potentials in the presence of two lighter quarks qq and 
provide parameterizations of these results by continu¬ 
ous functions. In section EYi we use these parameteri¬ 
zations in model calculations and check for the existence 
of bound states, which would indicate the existence of 
tetraquarks. We conclude in section [V| 


II. MODELING THE bb INTERACTION IN THE 
PRESENCE OF TWO LIGHT QUARKS qq 


while for a quark-quark pair q^^Aqj 

C = -^Tr(A“AA“'^Al') (3) 

a 

with the Gell-Mann matrices A“, a = 1,..., 8. For exam¬ 
ple, Kab = Sab /v^ describes the qq color singlet, while 
Aab = iAB 3 /is one of three independent possibilities 
to realize a qq color triplet. In Table |lj the resulting val¬ 
ues for C for the qq singlet and octet and the qq triplet 
and sextet color orientations are listed. 


color 

qq singlet 

qq octet 

qq triplet 

qq sextet 

orientation 

1 

8 

3 (and 3) 

6 (and 6) 

C 

+ 16/3 
(attractive) 

-2/3 

(repulsive) 

+8/3 

(attractive) 

-4/3 

(repulsive) 


In the following, we discuss quark model expectations 
regarding the qualitative behavior of a qqbb four-quark 
systems, where q denotes either a light u, d^s or c quark 
[49| . In particular, we are interested in the bb interaction 
in the presence of two light quarks qq. The qualitative 
expectations are confirmed by corresponding lattice QCD 
results, which are discussed in section [TTTl The main pur¬ 
pose of these model considerations is to motivate a suit¬ 
able fit function for the lattice QCD bb potential results, 
which is used in section IV in the Schrodinger equation to 
check whether and in which channels bound four-quark 
states, i.e. tetraquarks, exist. 


A. The quark-antiquark / quark-quark potential at 
small separations 


In the original quark model |23j . the quark- 
antiquark and the quark-quark (or equivalently anti- 
quark-antiquark) potentials at small separations r = 
\ri — Tjl are dominated by one-gluon exchange similar 
to the Fermi-Breit interaction, 
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16 si • Sj 
Smiirij 



Table I: The color factors C for the qq singlet and octet and 
the qq triplet and sextet color orientations. 

Lattice QCD confirms that the static color singlet po¬ 
tential at small separations r can be described reason¬ 
ably well by one-gluon-exchange (cf. e.g. [26l[27], where 
a matching of lattice QCD and perturbative results is 
done). At larger separations, it becomes linear with cer¬ 
tain l/r-corrections due to string vibrations |28] . One 
can crudely estimate as appearing in Q by consider¬ 
ing the color singlet qq. In that case C = -1-16/3, while 
string vibrations lead to « —'KjVlr at intermediate 
separations, resulting in as ~ 7r/16. While this estimate 
is most appropriate for static quarks, as is expected to 
be somewhat larger for quarks of finite mass PH US] . 

The only spin dependent term in 0 is the hyperfine 
interaction proportional to • Sj, which is pathological 
in the original quark model due to the Dirac delta (cf. 
eq. 0). In the relativistic quark model, however, this 
interaction is smoother and, hence, well behaved P^I25j . 
Clearly, the interaction is weaker for a spin triplet than 
for a spin singlet. 

To summarize, whether the potential between a quark 
and another quark or antiquark is attractive or repulsive 
depends on their color orientation. For small separations 
it is approximately Coulomb-like with the color factors 
C collected in Table |lj The hyperfine term enhances the 
interaction for a spin singlet and decreases it for a spin 
triplet. 


{i,j are the (anti)quark indices, r^, and denote 
their positions, spins and masses, respectively). Since 
we are exclusively interested in ground states, we have 
specialized eq. Q to angular momentum 1 = 0. The 
quark model has been improved (cf. e.g. pU [25]), but 
maintains its main ingredients. C depends on the color 
orientation of the (anti)quarks, which can be specified by 
a 3 X 3 matrix A. For a quark-antiquark pair qiAqj 

C = -f ^Tr(A“AA“Al'), (2) 


B. Qualitative discussion of the qqbb system 

For the particular case of the qqbb system, where the bb 
pair is significantly heavier than the light qq pair, we uti¬ 
lize the Born-Oppenheimer approximation |16j : for the 
light quarks, the heavy antiquarks can be regarded as 
static color charges; once the energy of the light quarks 
is determined, it can be used as an effective potential for 
the heavy antiquarks. We assume that at small bb sep¬ 
arations r, the b quarks interact according to the quark 
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model discussed in section IIA while at larger separa¬ 
tions their interaction is screened by the light quarks, i.e. 
the four quarks form two rather weakly interacting B) 


(s,c) 


mesons denotes either a B, B*, Bg, B*, Be or B* 

meson). 


Expectations for the bb interaction at small separations r 



Figure 1: (Color online) At small bb separations r, the heavy 
antiquarks bb form an antidiquark, which corresponds to a 
color triplet. There is essentially no screening of the bb inter¬ 
action due to the much farther separated light quarks qq. 


• The spin interaction of the b quarks is quite small 
and can possibly be ne^cted, since it is propor¬ 
tional to llrrib'^ (cf. eq. Q). 

• In case of a bound qqbb state, i.e. a tetraquark, the 
antiquarks bb are expected to be in a color triplet 
3, which is attractive, and not in a color sextet 6, 
which is repulsive (cf. also Table |T|. In other words, 
at small separations r, the antiquarks bb form an 
antidiquark as depicted in Figure [l] 

• Because the complete four quark system qqbb nec¬ 
essarily forms a color singlet, the light quarks qq 
must be in a color antitriplet 3. 

• Since this color antitriplet is antisymmetric, and 
since the light quarks qq are assumed to be in a spa¬ 
tially symmetric s-wave, the Pauli principle implies 
a symmetric spin-flavor structure. This can either 
be a spin singlet with an antisymmetric flavor com¬ 
bination or a spin triplet with a symmetric flavor 
combination. Indeed, when studying light u and d 
quarks in the presence of two static antiquarks us¬ 
ing lattice QCD, two attractive channels have been 
found [niMiisn]. As expected, these are a (spin) 
scalar isosinglet (j = 0, / = 0, where j denotes 


the spin of the light quarks qq) and a (spin) vec¬ 
tor isotriplet {j = 1, I = 1). The scalar isosinglet 
is more attractive, as expected from the hyperfine 
interaction in eq. (§, i.e. the lattice QCD results 
confirm the qualitative quark model expectations. 

• When studying two identical light quarks qq = ss 
or qq = cc, which are symmetric in flavor, the only 
attractive channel is a spin triplet. However, it is 
conceptually interesting to consider two hypothet¬ 
ical degenerate flavors with the mass of strange or 
charm quarks and then also investigate spin singlets 
with flavor structure qq = — s^^^s^^^)/-\/2 

and qq = — c^‘^^c^^^)/V2. 


Expectations for the bb interaction at large separations r 

• At large separations r, screening of the bb interac¬ 
tion is expected due to the light quarks qq, as illus¬ 
trated in Figure^ When the bb separation is larger 
than around two times the radius of a B(^s c) m^son, 
there is essentially no overlap between the wave 
functions of the light quarks and, consequently, the 
bb interaction practically vanishes. 

• The more massive the light quarks are, the more 
compact their wave functions in the c) mesons, 
as shown in Figure [^a), (b) and (c) and, thus, the 
stronger the screening. In other words, the corre¬ 
sponding bb potential becomes more and more nar¬ 
row and will at some point not anymore be able to 
host a bound state. Consequently, for a sufficiently 
heavy pair of light quarks qq the screening should 
prevent the formation of qqbb tetraquarks. 


Quantum numbers of possibly existing qqbb tetraquarks 

We study exclusively states which correspond for large 
bb separations to pairs of mesons in a spatially 

symmetric s-wave. Therefore, the parity of these states 
is positive, i.e. P = -\- (the product of the parity quantum 
numbers of the two mesons, which are both negative). 

As argued above, the two antiquarks bb are expected 
to be in an antisymmetric color triplet. Since their flavor 
is symmetric, their spin jh must also be symmetric due to 
the Pauli principle, i.e. jt = 1. Similarly, for an antisym¬ 
metric qq flavor combination, i.e. qq = (ud — du) /v^, 
j — 0, while for symmetric flavor combinations, i.e. 
qq € {uu, {ud -|- du)/y/2,dd , ss , cc}, j = 1. The total 
spin J of the qqbb system is the combination of j and jh- 

Altogether, the possibly existing qqbb tetraquarks we 
are going to investigate have the following quantum num¬ 
bers: 

• qq = {ud — du) / 

/(J^) = 0(1+). 
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i X 

(d) up / down 

light 

quarks 



(b) 



Strange 

light 

quarks 



(C) 


i X 

charm 

light 

quarks 





from the QCD Lagrangian (cf. 


and section III). Of 


course, these potentials can be obtained only for a limited 
number of discrete separations r. Therefore, a suitable 
fit function is required, to interpolate between the lat¬ 
tice QCD results and also to extrapolate beyond them. 
This fit function is based on the qualitative expectations 
discussed above and will be used in the Schrodinger equa¬ 
tion in section m where we determine whether and in 
which channels bound four-quark states exist. 

For two heavy antiquarks bb inside a cloud of two 
light quarks qq, i.e. at small bb separations, we expect a 
Coulomb-like potential of order —2as/3r « —7r/24r cor¬ 
responding to a color triplet. At larger separations r, the 
potential will be screened by the light quarks qq. This is 
due to the decrease of the wave function ip of each of the 
light quarks with respect to their separations from the 
heavy antiquarks. One expects this decrease to follow an 
exponential of a power of r, i.e. ip oc exp(—(r/d)^), where 
d roughly describes the size of each of the bq systems, i.e. 
the size of a ^(Vc) niGSon <0.5fm. The parameter p 
characterizes the radial profile of the light quark wave 
function inside the ^(Vc) ™6Son. Assuming the qb inter- 


action inside the B; , meson is dominated by a linear 
confining potential, one can estimate the parameter p. In 
the case, where the quark q is rather heavy, e.g. q = c, 
the corresponding non-relativistic Schrodinger equation 
is solved by Airy functions, resulting in p = 3/2. A sim¬ 
ilar but relativistic treatment for a lighter quark yields 
p = 2 instead. 

These considerations suggest the following fit function 
for lattice QCD bb potential results: 


V{r) = -“exp(-(^)'’) +Fo, (4) 


r y 


Figure 2: (Color online). At large bb separations r, the qqbb 
system is essentially a system of two mesons. The 

color charge of each of the antiquarks b is almost completely 
screened by one of the light quarks q. (a) q £ {w, d}. (b) q = s. 
(c) q = c. 

• qq £ {uu, {ud + du)j'/2, dd}: 
/(J^)e{i(0+),i(i+),i(2+)}. 

• qq& {ss, cc}: 

/(j^)e{o(o+),o(i+),o(2+)}. 

C. Fit function for lattice QCD bb potential results 

Using lattice QCD, one can compute bb potentials in 
the static limit (i.e. for rri}, —>■ cxd) from first principles, i.e. 


where it is expected that a « 2as/3 « 7r/24 « 0.13, 
(i<0.5fm and p « 1.5... 2.0. The constant Vb is nec¬ 
essary to account for twice the mass of the static-light 
meson. As will be demonstrated in the following section, 
this fit function is consistent with lattice QCD results 
and the crude quantitative expectations for a, d and p 
are fulfilled. 


III. LATTICE QCD COMPUTATION OF THE bb 
INTERACTION IN THE PRESENCE OF TWO 
LIGHT QUARKS qq 

To determine the effective bb potential quantitatively, 
we use lattice QCD and consider the limit of infinitely 
heavy b quarks, i.e. the static limit. The first lattice 
computations of such potentials have been performed in 
the quenched approximation (cf. e.g. |5TH55] ). Recently, 
also computations with dynamical sea quarks have been 
performed [l3llIl[29ll3Ql|36l[3a. In this work, we ex¬ 
tend our previous computations for light quark combi¬ 
nations qq € {{ud — du)j^/2 , uu,{ud + du)/'/2,dd} 
[291130] by similar computations with strange and charm 
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quarks, i.e. 99 € , ss , — 

c( 2 )c(i))/-y /2 , cc}. 


A. Lattice QCD setup 


We have performed computations using two ensembles 
of gauge link configurations generated by the European 
Twisted Mass Collaboration (ETMC) with 2 dynamical 
quark flavors. The quark action is Wilson twisted mass 
tuned to maximal twist, while the gluon action is tree- 
level Symanzik improved. Most importantly, this guar¬ 
antees automatic 0{a) improvement of spectral quanti¬ 
ties, i.e. discretization errors in the resulting bb potentials 
appear only quadratically in the lattice spacing a. Infor¬ 
mation about these ensembles is collected in Table El 
Further details, in particular regarding their generation, 
can be found in [3H1135] . 


p 

size 

fJ-i 

a in fm 

niT, in MeV 

configurations 

3.90 

24® X 48 

0.0040 

0.079 

340 

480 

4.35 

32® X 64 

0.00175 

0.042 

352 

100 


Table II: Ensembles of gauge link configurations used for the 
computation of bb potentials (/3: inverse gauge coupling; size: 
number of lattice sites; /rj: bare u/d quark mass in lattice 
units; a: lattice spacing; m^: pion mass; configurations: num¬ 
ber of gauge link configurations used). 


Eor bb potentials in the presence of two light quarks qq 
with 9 S {u,d}, we reuse our lattice QCD results from 
P51I5U] . which were obtained using the ensemble with the 
coarser lattice spacing a « 0.079 fm. For q G {s, c}, the bb 
interaction is screened at significantly smaller bb separa¬ 
tions (cf. the discussion in section IIB and Figure]^. To 
be able to resolve the corresponding potentials properly, 
we decided to use for flavor combinations qq G — 

g(2)g(i))/y2 , ss , , cc} another 

ensemble with a finer lattice spacing a » 0.042 fm. Al¬ 
though the physical extent of the lattice for this ensemble 
is much smaller than for the other one, this should not 
introduce significant finite volume effects at the rather 
small separations we are interested in. 


Note that for both ensembles, the u/d quarks are un- 
physically heavy, corresponding to a pion mass « 
340 MeV. Moreover, there are no s and c sea quarks, 
i.e. our lattice QCD results are obtained in a partially 
quenched approximation. For the computation of bb po¬ 
tentials in the presence of light s and c quarks, we also 
use a much smaller number of gauge link configurations. 
The reason is that the propagators of the heavier s and 
c quarks introduce less statistical noise than those for 
lighter u/d quarks. 


B. Lattice QCD computation of bb potentials 


We determine bb potentials in the presence of two light 
quarks qq from the exponential decay of temporal corre¬ 
lation functions, 

C(t,|ri-r2|) = (Q|O^(t)O(0)|Q) (5) 

of four-quark creation operators 

0{t) = {CT)AB{Ct)cD 

(Qc(ri)9A^(ri)) (Qr>(r2)9B^(r2)) ( 6 ) 

at sufficiently large tmin ^ t < tmax- Here Q denotes a 
static antiquark operator approximating a b quark, 9 is a 
light quark operator, A,B,C,D are spin indices, (I), (2) 
are flavor indices and C = 7072 is the charge conjugation 
matrix. For the static antiquarks, the only relevant vari¬ 
able is their separation. Their spin components can be 
combined with f e_{(l - 70)75, (1 - 7 o) 7 j}, J = 1,2,3, 
where the resulting bb potential does not depend on which 
F matrix is chosen. The spin components of the two light 
quarks can be coupled in 16 independent ways via T, 
which should be an appropriately chosen combination of 
7 matrices to realize definite quantum numbers \jz\ (an¬ 
gular momentum with respect to the axis of separation), 
P (parity) and Px (behavior under reflections across an 
axis perpendicular to the axis of separation). For a more 
detailed discussion of symmetries and quantum numbers, 

cf. m- 

Note that the creation operators Q, when applied to 
the vacuum |Q), do not only generate definite quantum 
numbers {\jz\,P, Px), but also a structure resembling two 
mesons separated by r = |ri — r 2 |. Such operators 
should be well suited to excite the ground state of the 
corresponding (\jz\, P, Px) sector, in particular for large 
QQ separations r, where one expects two weakly inter - 


acting mesons (cf. the discussion in section 

Note, however, that the arrangement of the four quarks 
qqQQ in the ground state is decided by QCD dynamics, 
i.e. automatically realized in the lattice result accord¬ 
ing to QCD and not by the structure of the employed 
creation operators. For example, in recent lattice QCD 
work on tetraquark candidates, it has been demonstrated 
that operators similar to generate significant overlap 
to a variety of different four-quark structures, including 
mesonic molecules, diquark-antidiquark pairs and two es¬ 
sentially non-interacting mesons gnmi]. 

In previous computations [131 ESI EQ] , we have consid¬ 
ered light quarks 9 G {u, d} (due to technical reasons, 
the quark mass mu,d was chosen unphysically heavy cor¬ 
responding to a pion mass « 340 MeV; cf. also the 
first line in Table[ll| . We studied the scalar isosinglet with 
antisymmetric spin j = 0 and flavor 99 = {ud — du)/^/2 
(in the following denoted as the scalar u/d channel), as 
well as the vector isotriplet with symmetric spin 9 = 1 
and flavor qq G {uu, {ud + du)/ v^, dd} (in the following 


IIBl. 
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denoted as the vector u/d channel), which are the two 
attractive channels between ground state mesons {B and 
B*). Note that the scalar u/d channel was found to be 
more attractive than the vector u/d channel, as expected 
from quark model cons iderations (cf. eq. Q and the dis¬ 
cussion in section IIB). 

In this work, we extend these computations to heavier 
pairs of light quarks qq = ss and qq = cc. For these sym¬ 
metric flavor combinations, the only attractive channel 
for two ground state mesons {Bi/}) is the vector chan¬ 
nel, i.e. with light quark spin j = 1. It corresponds to 
r = (I -I- 7 o) 7 j, j = 1,2, 3, in the creation operator §. 

To be able to study also the scalar channel, i.e. j = 0, 
with strange and charm quarks, we consider two hypo¬ 
thetical degenerate flavors with the mass of the strange 
or the charm quark, which allow to form antisymmet¬ 
ric flavor combinations qq = 

and qq = It corresponds to 

r = (I -I- 70)75 in the creation operator Q. 

For further details regarding the lattice QCD compu¬ 
tation of bb potentials, we refer to [29l |30]. Examples 
for qq = {ud — du)/\/2 (scalar u/d channel) and for 
qq € {uu, (ud + du)/\/2, ddj (vector u/d channel) are 
shown in [13], Figure 1. 


C. Fitting eq. ([^ to lattice QCD bb potential 
results 


To describe the lattice QCD bb potential results I/^®'*(r) 
by continuous functions, we perform uncorrelated 
minimizing fits of eq. Q, i.e. we minimize 


= 


E 


/ F(r)-F‘^t(r) 


(7) 


with respect to the parameters a, d and Vp, w hile keeping 
p = 2 fixed (cf. the discussion in section IICI [50|. 
denote the corresponding statistical errors. 

We perform these fits for the scalar u/d, the vector u/d, 
the scalar s, the vector s and the scalar c channel. The 
lattice QCD bb potential of the remaining vector c chan¬ 
nel is, however, strongly screened and consistent with 
yiat^r^ = 0 for r > 2a. Such results are not sufficient to 
perform a stable fit. 

To investigate and quantify systematic errors, we do 
not only perform a single fit for each of the mentioned 
five channels, but a large number of fits, where we vary 
the following parameters: 


• The range of temporal separations tmin < ^ < ^max 
of the correlation function C(t, r) (eq. ([^) at which 
yiat(j,) jg according to: 


^max ^min ^ ^7 

— for u/d channels: 

4(2 ^ tmin; tmax ^ 9^; 


— for s and c channels: 

10 a 

E ^min E 14a, t-max ^ 19a 

(small tmin might lead to a contamination by ex¬ 
cited states; large tmin and tmax drastically increase 
statistical errors). 

• The range of spatial bb separations < r < r„iax 
considered in the minimizing fit Q, according 
to: 


— for the vector u/d channel: 

^min — 2a [51] : 

— for all other channels: 

^min € '[2a, 3a}, 

— for u/d channels: 

^max € [^a, 9a, 10a}, 

— for s and c channels: 

^max € {7a, 8a} 

(yiat(j,) small r < 2a are expected to suffer from 
sizable lattice discretization errors, while F*^*(r) at 
large r is essentially a constant, i.e. has little effect 
on the relevant fit parameters a and d). 


For each of the fitting parameters a, d and Vq, we con¬ 
struct a distribution by considering the results of all the 
above listed fits weighted by exp(—y^/dof) with from 
eq. (^. The central values of a, d and Vp are then de¬ 
fined as the medians of the corresponding distributions 
and the lower/upper systematic uncertainties are given 
by the difference of the 16th/84th percentiles to the me¬ 
dians (in the case of a Gaussian distribution, an uncer¬ 
tainty defined in this way would correspond to its width, 
i.e. la). Since in general the distributions are asymmet¬ 
ric, the systematic uncertainties are asymmetric as well. 
For more details regarding this method of estimating sys¬ 
tematic errors we refer to [42] . 

Finally, to include statistical errors, we compute the 
jackknife errors of the medians oi a, d and Vp and add 
them in quadrature to the corresponding systematic un¬ 
certainties. 

To illustrate this error estimation procedure, we show 
in Figure example histograms representing the distri¬ 
bution of a and d for the scalar u/d channel. The green, 
red and blue bars correspond to the systematic, statisti¬ 
cal and combined errors, respectively. In the following, 
we will always use and quote the combined errors repre¬ 
sented by the blue bars. 

The hnal results for a and d are collected in Table m 
Note that within errors they agree with the model con¬ 
siderations and crude quantitative expectations discussed 
in section jn] We do not list results for Vp, since it is an 
irrelevant constant corresponding to twice the mass of a 
static-light meson. The fit function Q with the parame¬ 


ter sets from Table ]III 
are shown in Figure [4 
qualitative expectations discussed in section ]ITI3 


and the corresponding error bands 
Clearly, these results confirm the 
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a 



Figure 3: (Color online). Histograms used to estimate sys¬ 
tematic errors for a and d for the scalar ujd channel (green, 
red and blue bars represent systematic, statistical and com¬ 
bined errors, respectively). 


(1) The screening of the bb interaction is stronger for 
heavier light quarks qq. 

(2) The scalar channels are more attractive than the 
corresponding vector channels. 



Spin 

a 

d in fm 

(ud 

— du )/ 

scalar 

Q or+0.04 
U.OO_0.04 

U.+^_0.08 

uu, {ud + du)/\/2, dd 

vector 

u.zy o.o6 




scalar 

r\ Qy+0.08 

' -0.05 

0.201°:;° 


ss 

vector 

'^•-*-°-0.02 

0 18+°'^^ 
'^•-*-°-0.05 


-c(")cW)/V2 

scalar 

0 10+0.12 
'^•-‘-^-0.07 

0 


Table III: Parameters a and d obtained from minimizing 
fits of Q to lattice QCD bb potential results. 


IV^ DEPENDENCE OF THE EXISTENCE OF 
qqbb TETRAQUARK STATES ON THE LIGHT 
QUARK MASS 


In [13], we have found evidence for a bound state in 
the scalar u/d channel, i.e. the existence of a qqbb = 
udbb tetraquark. For heavier quarks qq, the effective bb 
potentials are less attractive. This has qualitatively been 


anticipated in section |TT] and quantified in section III (in 
particular cf. the resulting values for a and d in Table [TTl| 
and the plots in Figure]^. Thus, for a sufficiently heavy 
pair of light quarks qq we expect that the qqbb system 
will not anymore be able to form a bound state. In the 
following, we investigate whether this is already the case 
for strange and/or charm quark masses. We also study 
the vector channels. 


A. The bb Hamiltonian 


We define U{r) = V{r)\vg=o,p =2 with V{r) from of eq. 
Q. U{r) with a set of fit parameters a and d from Ta¬ 


ble III corresponds to the ground state energy of a qqbb 


4-quark system in a specific channel minus the energy of 
a pair of far separated mesons. Thus, the corre¬ 

sponding Hamiltonian for the relative coordinate of the 
bb quarks is 


H = ^+2m^^ + C/(r), (8) 

where = itihI^ is the reduced mass. At large sep- 
arations, each b quark carries the mass of a me¬ 
son because of screening, and thus ttih = . At 

(s,c) 

small separations, mn = mb could be more appropriate. 
Throughout this section, we always consider two choices, 
mH = mB^^ and mH = mb, which yield qualitatively 
identical results. Note that any dependence on the heavy 
b spins is neglected, because V (r) has been computed in 
the static limit ^ oo. Since the b quarks are quite 
heavy, we expect the static limit to be a reasonable ap¬ 
proximation. 

In classical mechanics, the bb separation r would vanish 
for the ground state, but after quantizing the system, a 
bound 4-quark state {E < 2mH) may not exist anymore. 


B. An analytical estimate for qqbb binding 


In |I3j , we have derived an approximate analytical rule 
for the existence/non-existence of a bound qqbb state us¬ 
ing the Bohr-Sommerfeld quantization condition. If 


> 

- 128 X 2i/p(r(I-bl/2p))2 


(9) 


is fulfilled, there should be at least one bound state. The 
right hand side of this rule has a rather moderate depen¬ 
dence on the exponent p. For example, when p increases 
from the expected values of 1.5 to 2.0 (cf. section IICl, 
the right hand side only changes from 0.55 to 0.60. Thus, 
the existence of a bound state mainly depends on the 
product of parameters pad. 

With the medians for the parameters a and d (cf. Ta¬ 
ble III), we determine the left hand side of eq. ((^. For 
the reduced mass, we use both uih = ms^^ {ms = 












































Figure 4: (Color online), bb potentials in the presence of two lighter quarks qq {qq flavor: up/down in green, strange in blue, 
charm in red; qq spin: j — 0, i.e. scalar, in the upper line, j = 1, i.e. vector, in the lower line). The plotted curves with the error 
bands correspond to eq. (fill with the parameter sets from Table III Vertical lines indicate lattice separations r = 2a, 3a,... of 


lattice QCD potential results V‘“(r) used to generate the parameter sets from Table III via x minimizing fits. 


5279 MeV, tub, = 5367 MeV, tob, = 6276 MeV gS]), 
which is certainly a good choice for large bb separations, 
and ^ = mb/2 {mb = 4977MeV, from quark models gi]). 
which might be more appropriate for small bb separations 
(cf. the discussion in section IV AI. The results for fiad 
for the u/d, s and scalar c and vector channels are col¬ 
lected in Table IV For the scalar u/d channel, there is 
strong indication for the existence of a tetraquark (i.e. 
Iiad ^ 0.60), which confirms our findings from [13) . For 
the vector s channel and for charm quarks, bound qqbb 
states are not expected (i.e. ^ad ^ 0.60). For the vector 
u/d and the scalar s channel the situation is less clear. A 
more rigorous and quantitative analysis is needed, which 
is part of the following section. 


C. Numerical solution of the Schrodinger equation 


fiad 



Spin 

mu = 

mn = rrib 

{ud — du)/-\/2 

scalar 

1.97 

1.86 

uu, {uddu)/dd 

vector 

0.60 

0.57 

( s ( 1 ) s ( 2 )_ 5 ( 2 ) s ( 1 ))/^ 

scalar 

0.74 

0.69 

s$ 

vector 

0.44 

0.41 

(c(l)c(2) -c(2)c(l))/y2 

scalar 

0.34 

0.27 


Table IV: Values for fictd, which represent the left hand side 
of eq. @. Values >0.60 (right hand side of eq. § for p = 2) 
point towards the existence of a bound qqbb state, while values 
< 0.60 are an indication against the existence of such a state. 


expected in an s-wave, for which the radial equation is 

= {E-2mH)Rir) ( 10 ) 

=Eb 


To investigate the existence of a bound state more 
rigorously, we solve the Schrodinger equation with the 
Hamiltonian ([^ numerically. The strongest binding is 


with the wave function ip = ip{r) = R{r)/r. If Eb = E — 
2mH < 0, —Eb can be interpreted as the binding energy. 
We proceed as explained in m and solve this equation by 
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imposing Dirichlet boundary conditions R{r = oo) = 0 
and using 4th order Runge-Kutta shooting. 

For the scalar u/d channel, the lowest eigenvalue Eg < 
0 , which implies the existence of a bound four-quark 
state. For all other channels, i.e. the vector u/d and 
the s and c channels. Eg > 0 , i.e. the corresponding qqbb 
tetraquarks will most likely not exist in these channels 
[5^ . These hndings conhrm the analytical estimates ob¬ 
tained in the previous subsection (eq. Q and Table [lV| . 

The central value and the combined systematic and 
statistical error for the binding energy Eb of the 
tetraquark state in the scalar u/d channel is obtained 
by the method discussed in section III C (gener ating a 
distribution for Eb from the fits listed in section III C): 


Eb = —90^42 MeV (for ttih = tub), (H) 

Eb = -93tg MeV (for uih = mb). (12) 


These binding energies are roughly twice as large as 
their combined systematic and statistical errors. In 
other words, the confidence level for this udbb tetraquark 
state is around 2a. The corresponding histogram for 
mn = wb is shown in Figure 



Eg in MeV 

Figure 5: (Color online). Histogram used to estimate the 
systematic error for the binding energy Eb for the scalar u/d 
channel and mu = tub (green, red and blue bars represent 
systematic, statistical and combined errors, respectively). 



spin 

mn = msy y 

rriH = mb 

(ud — du)/-\/2 

scalar 

0.46 

0.49 

uu, {ud + du)/y/2, dd 

vector 

1.49 

1.57 

(s(l)5(2)_5(2)s(l))/y2 

scalar 

1.20 

1.29 

s$ 

vector 

2.01 

2.18 

(c(1)c(2)c(2)cW)/V2 

scalar 

2.57 

3.24 


Table V: Factors by which the mass rriH has to be multiplied 
to obtain a tiny but negative energy Eb . Factors <C 1 indicate 
strongly bound states, while for values 3> 1 bound states are 
essentially excluded. 


constant are plotted in the a-d-plane starting at a tiny 
energy Eb = —0.1 MeV up to rather strong binding, 
Eb = —100 MeV (gray dashed lines have been computed 
with mH = mBf^ gray solid lines with rriH = mb). 
The three plots correspond to u/d,s and c light quarks 
qq, respectively. Each fit of eq. (M to lattice QCD bb 
potential results (cf. the detailed discussion a bout s ys- 
tematic error estimation for a and d in section III CI is 


represented by a dot (red: scalar channels; green: vector 
channels; crosses: r^nin = 2a; boxes: = 3a). The 

extensions of these point clouds represent the systematic 
uncertainties with respect to a and d. If a point cloud is 
localized above or left of the isoline with Eb = —0.1 MeV 
(approximately the binding threshold), the correspond¬ 
ing four quarks qqbb will not form a bound state. A 
localization below or right of that isoline is a strong in¬ 
dication for the existence of a tetraquark. In case the 
point cloud is intersected by that isoline, the estimated 
systematic error is too large to make a definite statement 
regarding the existence or non-existence of a bound four- 
quark state. The big red and green bars in horizontal and 
vertical direction represent the combined systematic and 
statistical errors of a and d, as quoted in Table |III[ One 
can observe and conclude the following from Figure 

• There is clear evidence for a tetraquark state in the 
scalar u/d channel. 


To crudely quantify also the non-existence of bound 
four-quark states in the remaining channels, we deter¬ 
mine numerically by which factors the heavy masses tob 
in the Schrodinger equation (10) have to be increased to 
obtain bound states, i.e. tiny but negative energies Eb 
(the potentials U{r) are kept unchan ged, i.e. we stick to 
the medians for a and d from Table III I. The resulting 
factors are collected in Table [V] While the scalar s chan¬ 
nel is quite close to be able to host a bound state, the 
scalar c channel and the vector channels are rather far 
away, since they would require b quarks approximately 
1.6... 3.3 times as heavy as they are in nature. Note 
that the factors listed in Table El could also be rele¬ 
vant for quark models aiming at studying the binding 
of tetraquarks quantitatively. 

In Figure we present our results in an alterna¬ 
tive graphical way. Binding energy isolines EB(pt^d) = 


• The scalar s channel is close to binding/unbinding. 
A definite statement with our currently available 
lattice QCD data is not possible. 

• the scalar c and all vector channels do not host a 
bound four-quark state. 

These findings are consistent with the results presented 
above in Table US and Table |Vl 


V. CONCLUSIONS AND OUTLOOK 

In a previous publication [13] , we have found indication 
for the existence of a qqbb tetraquark with qq = (ud — 
du)/'/2 (i.e. in the scalar u/d channel). In this work, 
we have extended these studies by considering for qq not 
only w/d, but also heavier s and c quarks. In contrast to 
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qq = (ud-du)/\/2 qq= uu, (ud+du)/V2, dd 



qq = qq = ss 



qq = (c(1)c(2>-c(1)c(2))/^2 



Figure 6: (Color online). Binding energy isolines i5s(a,d) = 
constant in the a-d-plane for u/d, s and c light quarks qq, 
respectively (gray dashed lines: mn = e)' gray solid 

lines: mn = mb). The red and green dots represent the fits 
of eq. Q to lattice QCD bb potential results, while the red 
and green bars are the corresponding combined systematic 
and statistical errors. 

m, we have also investigated and quantified systematic 
uncertainties in detail. 

Our main results are the following: 

• We confirm the udbb tetraquark state in the scalar 
u/d channel predicted in our previous paper [13] 
with confidence level « 2a. The overall quantum 
numbers of this state are I{J^) = 0(1'''). 

• There is no bound four-quark state in the vec¬ 
tor u/d channel (/(J^) =€ {1(0+), 1(1'*'), 1(2+)}). 


Note, however, that we have been using unphysi- 
cally heavy u/d quarks (m^ « 340 MeV). Since de¬ 
creasing the light quark mass should enhance bind¬ 
ing, it will be interesting to explore in the future 
whether a bound four-quark state exists at physi¬ 
cally light u/d quark mass. 

• ssbb and ccbb tetraquarks, which correspond to the 
vector s and c channels (J^ S { 0 +,l+, 2 +}), do 
not exist. 

• It is of conceptual interest to introduce a hypo¬ 

thetical second s or c quark flavor. Then it is 
possible to also study the scalar s and c chan¬ 
nels, i.e. — s*^^^s*^^^)/-\/ 2)66 and — 

£;(2)g(i))/y^)55 systems (J+’ = 1+). While in the 
scalar c channel there is no bound four-quark state, 
the situation is less clear for s quarks. Improved 
lattice QCD results (less statistical errors, finer res¬ 
olution of bb separations) are needed before a defi¬ 
nite statement can be made. Binding in the hypo¬ 
thetical scalar s channel would indicate a fortiori 
binding for four-quark systems {{us — su)/V2)bb 
and {{ds — sd)/'/2)bb. Such light-strange chan¬ 
nels would then be highly relevant for experimental 
tetraquark searches. 

We consider these results to be important because they 
indicate both to experimental collaborations and to 
quark model phenomenologists which qqbb tetraquarks 
are expected to exist and which are not. 

To supply data for future quark model studies of 
tetraquarks, we also provide parameterizations of the po¬ 
tential of two static antiquarks bb in the presence of two 
lighter quarks qq, where qq € {{ud — du)/\/2 , uu, {ud + 
du)/V2,dd , — s^'^^s^^^)/\/2 , ss , — 

c(2)£,(i))/v^}. Moreover, we have determined quanti¬ 
tatively for these channels by which factor the heavy 
quark or meson mass rriH has to be increased to obtain 
a tetraquark state. 

It is also interesting to compare our findings to other 
groups studying the same or similar systems using, how¬ 
ever, different theoretical approaches. For instance in 
m, in the framework of QCD sum rules, binding for fla¬ 
vors equivalent to udbb, usbb and s.sbb has been found, 
and no binding for doubly charmed tetraquarks. How¬ 
ever, these bound systems have = 0 “ and J+ = 1 “ 
different from our results. Another example using the 
Dyson-Schwinger framework is [48j . where a tetraquark 
composed of four charm quarks, i.e. cccc, has recently 
been predicted with a mass significantly lighter than 
2m^^. In principle our static antiquarks can also be con¬ 
sidered as a crude approximation of cc. Since we do not 
find a bound state for qq = cc, there seems to be a qual¬ 
itative discrepancy to our results, which would be inter¬ 
esting to understand and to resolve. 

As an outlook, it would be interesting to decrease the 
light u/d quark mass to their physical value, since this 
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should increase the radius of a i? meson, reduce screen¬ 
ing and, therefore, lead to a larger binding energy. As 
mentioned above, a tetraquark could then also exist in 
the vector u/d channel. Additionally, lighter u/d quark 
masses may also allow the study of light meson exchange 
interactions between the two B mesons. Because sim¬ 
ulations and computations at lighter u/d quark masses 
are computationally very expensive, we leave them for a 
future publication. 

Since there is a bound state for qq = {ud— du) /and 
possibly even for qq = it will be 

very interesting to investigate usbh (or equivalently dsbh) 
systems. This will, however, require additional compu¬ 
tations and also the implementation of certain modifica¬ 
tions in our analysis procedure. We plan to study such 
flavor combinations in the near future. 

Another interesting, but very challenging task, is to 
include corrections due to the heavy bb spins. While in 
principle it is possible to compute such corrections us¬ 
ing lattice QCD (cf. [HJ 03], where this has been pi¬ 
oneered for the standard static quark-antiquark poten¬ 
tial), in practice we expect this to be extremely hard for 
qqbb systems. Therefore, a more promising and realis¬ 
tic approach seems to include such spin-dependent in¬ 
teractions in the Schrodinger equation, which will result 
in a coupled channel differential equation. We are cur¬ 


rently in the process of exploring this approach, where 
first promising qualitative results have recently been pre¬ 
sented at a conference 05] . 

Once these techniques are fully developed for qqbb sys¬ 
tems, it will be most interesting to extend them to qqbb 
systems and to study the crypto-exotic bb tetraquark can¬ 
didates observed by the BELLE collaboration [3] . 
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